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Abstract 

Strong interaction level shifts and widths in S~ atoms are analyzed by using 
a X nucleus optical potential constructed within the relativistic mean field 
approach. The analysis leads to potentials with a repulsive real part in the 
nuclear interior. The data are sufficient to establish the size of the isovector 
meson-hyperon coupling. Implications to £ hypernuclei are discussed. 
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I. INTRODUCTION 



The relativistic mean field (RMF) approximation [1] is an interesting alternative to more 
conventional approaches to the nuclear many-body problem. Numerous studies have demon- 
strated the usefulness of the RMF approach in describing the behavior of nucleons in the nu- 
clear medium. Recent RMF calculations of A- and multi A hypernuclei [2-4] have shown that 
it can be successfully extended to the more general baryon- nucleus systems [5-9] . Hyperon- 
nucleus scattering has been studied as well [10]. In these works, the meson-hyperon coupling 
constants were either fitted to known hypernuclear spectra or obtained from the parameters 
of a nucleon shell model by using the constituent quark model. 

Calculations of S-hypernuclei [11] have often been based on analyses of S~ atomic data 
[12,13] in terms of attractive and absorptive S-nucleus potentials, yielding in the nuclear 
interior depth values 



However, a considerably shallower potential, particularly for ReV^ t , is indicated by fitting 
to S-hypernuclear spectra from (K~ ,7r + ) reactions [14]. Furthermore, a very recent phe- 
nomenological analysis of level shifts and widths in S~ atoms by Batty et al. [15] suggests 
that ReV p t is attractive only at the nuclear surface, changing into a repulsive potential 
as the density increases in the nuclear interior. Although the magnitude of the repulsive 
component cannot be determined unambiguously by the atomic data, the small attractive 
pocket of such potentials does not provide sufficient binding to form S-hypernuclei. This 
conclusion has been independently arrived at by Harada [16], folding a soft-core repulsive 
plus long-range attractive SA" interaction into the nuclear density. 

In view of the above seemingly conflicting observations, it is topical to use the RMF 
approach directly for constructing the S-nucleus optical potential by fitting to S~ atomic 
data. The aims of the present work are: 



ReV O pt(0) « 25-30 MeV, 
kVjj(O) « 10-15 MeV. 



(1) 



(2) 
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i) to check whether or not the RMF approach can reproduce the atomic data; 

ii) to establish constraints imposed by fitting to the data on the scalar and vector meson 
couplings to a S hyperon; 

iii) to establish the size of the isovector-vector meson-hyperon coupling; 

iv) to investigate the sign of the S-nucleus potential in the nuclear interior and whether 
or not £ hypernuclear bound states are expected to exist. 

The paper is organized as follows. Sect. 2 presents updated fits of a phenomenological 
density dependent (DD) optical potential. In sect. 3 we introduce the RMF optical potential 
for E, discuss the constraints imposed by the atomic data on its form and present results. 
These RMF results are compared with results of other approaches in sect. 4 which also offers 
a discussion and summary of our conclusions. 

II. PHENOMENOLOGICAL DENSITY DEPENDENT OPTICAL POTENTIAL 

Before comparing the E-nucleus potentials based on the RMF approach to those derived 
using the phenomenological analysis, we took the intermediate step of replacing the "macro- 
scopic" nuclear densities (such as 2-parameter Fermi distributions) used in the previous 
phenomenological analysis [15] by more realistic "microscopic" densities. The motivation 
for this replacement is twofold: first, we wish to use densities that are more consistent 
with the resulting RMF densities, and second, it is desirable to use densities that are as 
appropriate as possible outside the nuclear surface, for discussing £~ atoms. 

In order to appreciate the latter point, we have repeated the "notch test" of ref. [15] in 
a modified form to see which range of radii has a significant effect on the atomic levels. The 
sensitivity of the data to the potential around the radius is obtained by multiplying the 
best-fit potential by a factor: 

f=l-dexp[-(— ^) 2 ] (3) 
a N 
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representing a "notch" in the potential around r = Rn spread approximately over ±aN (in 
fm), whose relative depth is d. By varying the depth d and observing the changes in x 2 , 
the sensitivity (defined as the values of d that cause y 2 to increase by one unit) can be 
determined. For ojv the value of 0.5 fm was chosen. By scanning over R^, the radial region 
where the fit to the data is sensitive to the potential is determined. Figure 1 shows the 
results obtained, using "microscopic" densities as described below, by scanning R^ in steps 
of a nuclear diffuseness ao, i.e., Rn = i?o + a oA, with R = l.lA 1 ^ 3 fm and a = 0.5 fm. It is 
seen that the radial region which is "sampled" by the data is outside of the nuclear surface. 
This updates the similar results shown in Fig. 3 of ref. [15] using "macroscopic" densities. 

The method chosen to generate "microscopic" nuclear densities was to fill in single par- 
ticle (SP) levels in Woods-Saxon potentials separately for protons and neutrons. The radius 
parameter of the potential was adjusted to reproduce the rms radius of the charge distribu- 
tion (after folding in the finite size of the proton). The binding energy of the least bound 
particle was set equal to the corresponding separation energy. For N > Z nuclei the rms 
radius of the neutron density distribution was chosen to be slightly larger than that of the 
protons. By using this method, rather realistic density distributions are obtained in the 
region of the exponential fall-off. These densities are less reliable in the nuclear interior, but 
the E~ atom data are not sensitive to this region as is clearly seen from Fig. 1. 

Several of the fits of ref. [15] have been repeated using these SP densities. The S~ 
nucleus potential is written as 



2^(r) = -4tt(1 + 

m 



bo + Boi W) ] 



p(r) + 



5p(r) (4) 



>(0)' 

where \i is the S-nucleus reduced mass, m is the mass of the nucleon and p(r) = p n (r) + p p (r) 
is the nuclear density distribution normalized to the number of nucleons A. The isovector 
density distribution is given by 5p(r) = p n (r) — p p (r). The parameters b and B are given in 
units of fm and are obtained from fits to the £~ atom data. 

For the density-independent t e ffp potential (B — 0, B\ — 0), the results obtained with 
the SP densities were very similar to those obtained before [15], namely, values of x 2 P er 



degree of freedom, x 2 /F, were about 2.4 and the isovector part of the potential was not 
determined by the data. With the introduction of density dependence into the potential, 
substantial improvement in the fit to the data is observed. The best-fit values of x 2 obtained 
with SP densities are significantly lower than those obtained with the schematic parame- 
terization of the nuclear densities [15]. Moreover, the real part of the isovector potential is 
now found to be reasonably well-determined, presumably because the SP densities represent 
more reliably the differences between protons and neutrons outside of the nuclear surface. 

Results of two DD fits are summarized in table 1. Potential I is the analog of potential D 
or A' of ref. [15], where the density dependence is introduced only into the real isoscalar part 
via the parameter ReB . It is assumed throughout that lmbi=— lmb so that £~ particles 
are absorbed only on protons. All four parameters varied in the fit are phenomenological 
and it is seen that they are reasonably well-determined. The DD exponent a was held fixed 
during the parameter search and its value was subsequently varied between 0.1 and 1.1 and 
fits repeated. Potential II is analogous to the potential of ref. [15] for which the linear 
density terms were held fixed at values suggested by One Boson Exchange (OBE) models: 
&o=l-2-H0.45 fm, &i=— 0.45— i0.45 fm. In this case, three adjustable parameters, namely, 
the complex B and ReB 1 (ImB^- lmB ), were obtained from fits to the data. In both 
cases the fits are excellent and in both cases the variation of \ 2 with a was weak. The values 
shown in table 1 are typical examples and a could not be determined uniquely in the range 
of 0.4 to 0.9. 

From the parameter values in table 1 it is seen that whereas the potentials are attractive 
at low densities, they become repulsive as the density approaches nuclear matter densities, a 
feature that had already been observed [15]. Figure 2 shows the real and imaginary potentials 
for Si, Ca and Pb near the nuclear surface. It is seen that there is a small attractive pocket 
near the surface and that the real potential becomes repulsive toward the interior. The 
imaginary potential is purely absorptive; similar results for ReV^ t are obtained if ImV^ 
is made to saturate at nuclear matter density [17]. It should be mentioned that with the 
limited data available and the very different quality of various data, the isovector part of 
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the potential is determined solely by the experimental results for E atomic Pb. 



III. CONSTRUCTION OF THE RMF OPTICAL POTENTIAL 

The RMF formalism describes baryons as Dirac spinors coupled to scalar (a) and vector 
(u, p) meson fields. The underlying RMF model used here is based on the Lagrangian 
density of the form: 



£ = £ N + £ 



+ C v = - S>j ( 9 -f + ^■ / ,4"(t ; j a- ) S A , 



where 





( * S 






fi \ 


s = 





















and 



A 



(5) 
(6) 
(7) 



(8) 



(9) 



The standard form of C N can be found elsewhere [1]. The Lagrangian density £ s includes 
interactions of a £ hyperon with the isoscalar (a, a;) and isovector (p) meson fields, as well as 
with the photon generating in leading order the charged £ Coulomb field due to the nuclear 
charge distribution. The Euler-Lagrange equations lead to a coupled system of equations of 
motion for both baryon and meson fields. For the £ particle, and retaining only the timelike 
component p — for the vector fields in the mean (meson) field approximation, the resulting 
Dirac equation acquires the form: 



[-i a • V + V(r) + (3(M + S(r)) + ^couiomb(r) ]^(r) = E^(r) , 
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(10) 



where 



S(r) = g a ve(r) and V(r) = g u ^u(r) + g p shpo(r) (11) 

are the scalar and vector potentials, respectively. In eq. (11), I 3 = +1,0,-1 for S + , E°, 
E~, respectively. More details about the equations of motion and their solution are given 
in ref. [18]. 

For each particular nucleus we constructed the Schrodinger equivalent (SE) E-nucleus 
potential from the scalar (attractive) and vector (repulsive) Dirac potentials. The SE po- 
tential Vse{j : ) was then used as a real part of the optical potential in the calculations of £~ 
atoms: 

h^.w,- w + ^ + £6^ . (12 > 

It is to be stressed here that the isovector part of the real potential comes out quite naturally 
in the RMF approach as one allows for coupling the S with the p field (see eq.(ll) for V(r)). 

The imaginary part of the S-nucleus optical potential describes the conversion of £~ to 
A via the reaction T,~p — > An. The expression for the imaginary part of the potential (and 
consequently the S-hypernuclear width) was derived by Gal and Dover [13] and Dabrowski 
and Rozynek [19]. The values of the well-depth, — ImV^(O), were estimated to be between 
11.6 and 14 MeV in nuclear matter. For finite nuclei, a reduction from these estimates should 
be expected. In the present work a phenomenological imaginary potential of the form 

lmV l(r) = tp p (r) (13) 

was used, where the proton density p p was calculated in the RMF model and t was taken as 
a parameter to be determined by fitting the atomic data. The use of p p in the construction 
of the imaginary potential reflects the fact that the £~ conversion takes place exclusively 
on protons. 

For Cn the linear (L) parametrization of Horowitz and Serot [18] was used as well as the 
nonlinear (NL) parametrization NL1 of Reinhard et al. [20]. For Cy, the hyperon couplings 
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were characterized as in previous works [2-10] via the coupling constant ratios ccj = j^, 
(i — a, cu, p). In the case of the vector meson coupling ratio three values, namely 
1/3, 2/3 and 1, were used. The choice of a u = 1/3 was inspired by a number of early 
hypernuclear calculations [2,21-24] in which the strength of a hyperon (mainly A) coupling 
used was between 1/3 and 0.4. On the other hand, the 2/3 ratio follows from the constituent 
quark model. Finally, for a u = 1, the equality g^s = Q^n is motivated by a recent QCD 
sum rule evaluation for E hyperons in nuclear matter [25]. 

For each of the above choices of a w the ratio a a together with t from eq. (13) were fitted 
to the experimental atomic shift and width in Si. This particular atom was singled out due 
to its relatively accurate shift and width data among N = Z (isoscalar) core nuclei. The 
isovector part of the optical potential was then included and the ratio a p was determined 
by fitting the shift and width in Pb while holding the isoscalar parametrization fixed. In 
such a way, we found for each RMF model, whether L or NL, three sets of parameters 
that give, by construction, an excellent fit for Si and Pb. Typical results are illustrated in 
Fig. 3 for E nucleus NL optical potentials in Si. When going from q u = 1/3 to 1, ReV^ 
changes from attraction to repulsion in the nuclear interior. In addition, the imaginary part 
ImV^ t becomes more absorptive. The same holds for the linear model which, for a given 
value of a^, predicts more attraction for ReV^ t and less absorption for ImV^ than the 
NL parameterization yields, as is shown for a w = 2/3 in Fig. 4 where E optical potentials 
for Pb are compared with each other. We observed that the p - E coupling is determined 
unambiguously by the fit to the Pb data, with a value a p ~ 2/3 in all cases. 

Having determined the isoscalar as well as isovector parameters of V^ t by fitting the Si 
and Pb data we constructed SE optical potentials for all nuclei for which E~ atomic data 
exist. Whereas for q u = 2/3 and 1 we obtained reasonable values of x 2 without any further 
adjustments, we failed to get below x 2 = 26.5 for a w = 1/3. The results for the linear model 
are presented in Table 2. It is to be noted that we did not aim strictly at the lowest x 2 value 
but merely to get x 2 ~ 23, which is a typical value in previous DD fits [15]. Consequently, 
only the value of x 2 f° r the potential L3 from Table 2 represents a true minimum for this 
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particular parametrization. From this point of view a comparison of x 2 values listed in 
Table 2 might even be too favorable to the L3 potential. 

Although the S~ atomic data are of a poor quality and, moreover, determine the shape 
of Vf vt only around the surface and outside the nucleus as demonstrated by the "notch" test 
(Fig. 1), these data nevertheless significantly constrain the RMF parametrization of the £ 
nucleus optical potential. Not only that the value of the p-S coupling ratio a p ~ 2/3 holds 
unambiguously for all the parametrizations used, but in addition, a weak cu-T, coupling is 
almost certainly ruled out by fitting to the S~ atom data. Consequently, the data from XT 
atoms imply, within the RMF approach, a E-nucleus optical potential with a repulsion in 
the nuclear interior and a shallow attractive pocket outside the nuclear surface. For some 
light nuclei (Si, S, Mg), in contrast to Fig. 4 for Pb, the linear model with a w = 2/3 predicts 
that the repulsion in the interior goes back into a second shallow attractive pocket (less than 
3.5 MeV) at the origin, but this obviously cannot be tested by using £~ atom data. 

IV. DISCUSSION AND SUMMARY 

In this work we have improved on the phenomenological analysis [15] of S~ atoms in 
terms of a E nucleus DD optical potential by using nuclear density distributions which 
account more realistically for the nuclear surface region and outside of it. The earlier work 
[15] established that the isoscalar componenet of V^ t changes in the nuclear surface region 
from attraction to repulsion as one penetrates the nuclear interior, and that the isovector 
component is undetermined by the data. While the isoscalar component, in the present 
work, remains qualitatively unchanged, the introduction of these more realistic densities 
has a pronounced effect on the extraction of the isovector component of V^ t , which for 
S _ is reliably determined to be repulsive. It should be recalled that this extraction hinges 
almost exclusively on the shift and width data for Pb [26] which are the most accurate 
among all available £~ atom data to date, not just the N > Z data subset. More precise 
measurements, particularly for other N > Z atoms, are highly desirable. 
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We have also confirmed the results of ref. [15] that V^ t is determined by the E~ atom 
data only outside the nucleus. Nevertheless, the novel feature of V^ t becoming repulsive is 
found to occur at a region of space outside the nuclear radius, still where V^ t is determined, 
even if not very accurately, by the E~ atom data. To have a glimpse into the nuclear interior, 
a theoretical model is necessary, and in the present work we have used RMF as a working 
hypothesis. Although the RMF approach was applied before [7, 8] , for studying the binding 
of E hypernuclei, these studies were not constrained by any sound phenomenology since 
no E hypernuclear bound state has ever been clearly established [11], with the exception 
perhaps of a J n = + ,I = 1/2 |He bound state [27] which is a too light system to be 
used as input to RMF calculations. Furthermore, recent searches for E hypernuclear peaks, 
some of which had been reported on the basis of limited statistics, have yielded negative 
result [28, 29]. This leaves E~ atom data as the sole source of any possible E hypernuclear 
phenomenology. 

The application of the RMF calculational scheme to E^ atom data in terms of three 
coupling-constant ratios (a u , a a , a p ) has been successful in showing that very good quality 
fits, reproducing these data, can be made. The larger a w is (in the range to 1), the better 
is the fit. Of these fits the best ones, in the range a w ~ 2/3 to 1, indeed produce SE 
V^ t with a volume repulsion in the nuclear interior. The — 1 fit is an excellent one, 
reaching almost as low a x 2 value as those produced by the completely phenomenological 
DD optical potential fits of table 1. The observation that using relatively high values of ct u is 
distinctly superior to using relatively low values of a u within the RMF E-nuclear dynamics 
is in contrast to the situation in applying the RMF approach to A hypernuclei [2, 3, 6]. 
There, as stressed in ref. [30], a wider range of values for a u is roughly equally acceptable, 
since for any given value of a w the constraint of a A-nucleus potential well depth of about 
28 MeV attraction [31] can be satisfied by deriving an appropriate value for a a . In the 
absence of E hypernuclear data, no analogous constraint is operative for RMF E-nuclear 
applications. The volume repulsion obtained in this work for the isoscalar component of V^ t 
in fact precludes binding for E hypernuclei, at least for those based on core nuclei with a 

9 



small neutron excess, (N — Z)/A <C I. 

It is gratifying to conclude that both approaches, the phenomenological one and the RMF 
approach, agree with each other, for a comparable degree of fit, in producing an isoscalar 
repulsion in the nuclear interior and, for £~, an added isovector repulsion. We stress that 
this isovector repulsion, with a p ~ 2/3, was derived independently of the values assumed by 
the isoscalar coupling constants ratios. We recall that the common assumption of a pure F 
coupling for the p meson in SU(3), or equivalently of a universal p coupling to isospin, gives 
a p — 1, a value which is rather different from that determined by fitting to the S~ atom 
data. However, it is well known that in the OBE approach which unlike the RMF approach 
does use SU(3) guidelines for connecting the S — — 1 sector to the S = sector, the nuclear 
isovector single-particle potential is generated largely by other meson contributions 

than that due to a direct p exchange [32]. In the RMF approach, the p field is the only 
agency through which or can be generated, so that the statement a p ~ 2/3 is to 
be construed just as implying that ~ (2/3)Vi N \ where the baryon-nucleus isovector 
potential contribution is defined by 

vi B %-f/A , t B =^ W , (14) 

and T is the nuclear isospin operator with T 3 = (Z — N)/2. Since the nuclear isovector 
potential close to zero energy is estimated [33] as ~ 120 MeV, our estimate for v} S \ 
based on fitting £~ atomic Pb, is ~ 80 MeV. This value may be compared with the 
value \/ 1 (S) ~ 60 MeV [20], or \/ 1 (S) ~ 55 MeV [32], both estimates following Model D of 
the Nijmegen group [34] in agreement with the estimate derived from the phenomenology 
of 12 C(tf- tt*) reactions [35]. 

As in ordinary nuclear physics, the E-nucleus isovector potential cancels partly (for 
charged S ± ) the E Coulomb potential due to the nuclear charge distribution. For S + , we 
have checked that the attractive symmetry-energy contribution due to generally does 
not overcome the repulsive isoscalar contribution plus the repulsive Coulomb energy, so that 
it is unlikely to bind a S + in nuclei. For E°, where no symmetry energy or Coulomb energy 
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contribute, binding is precluded by the volume repulsion of the isoscalar E nucleus potential. 
These considerations have to be modified in very light hypernuclei where symmetry energy 
dominates over the Coulomb interaction (plus other charge-dependent effects due mostly 
to the mass differences within the £ charge triplet), leading to a better description of £ 
hypernuclei in terms of isospin eigenstates rather than charged £ states [35]. The only S 
hypernuclear system for which this dominance has been established by a detailed calculation 
[36] is |He where the S-nucleus isovector potential, represented as a Lane term of the form 
(14), is responsible for the binding of the / = 1/2, + state. Without this Lane term, |He 
would be unbound. 

For S~, the attractive Coulomb potential due to the nuclear charge distribution gives 
rise to an infinite set of £~ atomic states. Some of these states, however, for high Z nuclear 
cores, may be called nuclear states inasmuch as the corresponding wavefunctions are localized 
within the nucleus or in its immediate surface region. The energies and wavefunctions of 
these "nuclear" states generated by Vc can be approximated by noting that for a uniform 
charge distribution of a total charge Ze and radius R, Vc(r < R) is a harmonic oscillator 
potential 

V c = -V + ^Muj 2 r 2 , (15) 



Vo = -Za-^- , hu = \ - i (16) 

2 R V 3 Mc 2 R ' v ; 

with a depth V that grows as A 2 / 3 and an essentially A-independent oscillator spacing fuv, 

hoo = 3.25 ± 0.03 MeV between 40 Ca and 208 Pb for R = [(5/3) < r 2 >} l l 2 in terms of the 

nuclear charge rms radius < r 2 > 1 ^ 2 . For example, R = 4.49 fm for 40 Ca and R — 7.10 fm 

for 208 Pb. Requiring the rms radius of the S~ harmonic oscillator wavefunction to be less 

than or equal to R, one gets a series of nuclear radius values R^ 

R = 3.87 fm , R l = 5.00 f m , R 2 = 5.92 f m , R 3 = 6.71 f m , R A = 7.42 fm , (17) 

such that for Rn < R, the states of the first major shells up to and including N are "nuclear" . 
For 208 Pb, this suggest that these E harmonic oscillator states belonging to iV = 0,1,2,3 
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are "nuclear". This is confirmed by a numerical evaluation using the precise Vc, yielding 
the following binding energies 

Is : 20.35 MeV (18) 

lp : 16.92 MeV (19) 

Id : 13.63 MeV , 2s : 13.83 MeV (20) 

If : 10.58 MeV , 2p : 11.11 MeV (21) 

Including V^ t (L with = 2/3) in the binding-energy calculation, we find that these levels 
are displaced by 2 to 7.5 MeV to lower binding energies and also acquire strong-interaction 
widths. The Is level is very broad, about 36 MeV, whereas the widths of the other levels are 
relatively small and vary between 1.2 MeV to 2.6 MeV. The wavefunctions of these other 
levels are "pushed" by V^ t outside the nucleus. These states may be considered a special case 
of "Coulomb assisted states", a concept put forward in ref. [36]. A reaction which would 
excite preferentially a single £ value for £~ in 208 Pb, could locate these "nuclear" states 
for £ = 1,2,3. Their location will give a valuable information on the nature (attractive 
or repulsive) of V^ t in this mass region. However, since these states are "pushed" outside 
the nucleus, the cross sections to excite them by a nuclear reaction are not expected to be 
sizable. The chances of establishing a meaningful £ hypernuclear spectroscopy, therefore, 
are not particularly encouraging at present. 
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TABLES 



Table 1 - Parameters of DD E nucleus potentials. All parameters (except a) are in units 
of fm. Underlined parameters were held fixed during the fits. 



potential 


Re&o 


Im6o 


Re6i 


Re^o 


lmB 


ReSi 


a 


x 2 


X 2 /F 


I 


2.1 


0.50 


-1.0 


-4.3 








05 


16.7 


0.9 




± 0.8 


± 0.12 


± 0.4 


± 2.7 












II 




0.45 


-0.45 


-2.4 


-0.18 


-1.9 


08 


17.8 


0.9 










± 0.4 


± 0.19 


± 1.4 
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Table 2 - Parameters of the RMF S~ nucleus optical potentials. The linear parametrization 
from ref. [18] was used for the nucleonic sector, aj = g^/giN (i = cr, w, p), and t is defined 
in eq.(13). 



potential 




a<7 




t [MeV fm 3 ] 


x 2 


LI 


1.0 


0.770 


2/3 


-400 


18.1 


L2 


2/3 


0.544 


2/3 


-300 


23.9 


L3 


1/3 


0.313 


0.65 


-265 


26.5 
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FIGURE CAPTIONS 



Fig.l: Total x 2 values as a function of the relative notch depth d (cf. eq.3) superimposed 
on the best-fit DD optical potential I from Table 1 with a = 0.5 (cf. eq.4), for several 
positions Rn (Rn = Ro + ^oA) expressed via steps A of a nuclear diffuseness, with a 
radius R = I.IA 1 / 3 fm and a diffuseness a = 0.5 fm. 

Fig. 2: ReV^ (solid lines) and ImV^ (dashed lines) as functions of r for the best-fit DD 
optical potential I (cf. eq.4 and Table 1) for Si (Fig. 2a), Ca (Fig. 2b) and Pb (Fig. 2c). 
Arrows indicate the position of the corresponding nuclear rms radius. 

Fig. 3: ReV D p t (solid lines) and ImV^ (dashed lines) as functions of r for the RMF £~ 
optical potential NL in Si. Potentials in Figs. 3a, 3b and 3c correspond to the vector 
meson coupling ratios a w = 1/3, 2/3 and 1, respectively. 

Fig. 4: Comparison of the linear (L) (Fig. 4a) and nonlinear (NL) (Fig. 4b) RMF £~ 
optical potentials in Pb for the isoscalar vector meson coupling ratio = 2/3. The 
isovector vector meson coupling ratio is a p = 2/3 in both cases. 
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